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for  predicting  wave  reflection  coefficients  for  beaches,  revetments,  and  break¬ 
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sources  covering  a  wide  range  of  conditions  for  both  monochromatic  and  irregular 
waves.  The  work  was  carried  out  under  the  coastal  processes  program  of  the  U.S. 
Army  Coastal  Engineering  Research  Center  (CERC) . 

This  report  was  prepared  by  William  N.  Seelig,  Hydraulic  Engineer,  and 
John  P.  Ahrens,  Oceanographer,  both  of  the  Coastal  Processes  and  Structures 
Branch,  under  the  general  supervision  of  Dr.  R.M.  Sorensen.  J.  McTainany, 

Coastal  Oceanography  Branch,  provided  the  nonlinear  regression  analysis  used 
to  determine  empirical  coefficients  developed  in  this  report. 

Comments  on  this  publication  are  invited. 


Approved  for  publication  in  accordance  with  Public  Law  166,  79th  Congress, 
approved  31  July  1945,  as  supplemented  by  Public  Law  172,  88th  Congress, 
approved  7  November  1963. 


Colonel,  Corps  of  Engineers 
Commander  and  Director 
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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 


U.S.  customary  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x 

I0~3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angel) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins 

To  obtain  Celsius  (C) 

temperature 

readings 

from  Fahrenheit  (F)  readings. 

use  formula:  C  =  (5/9) 

(F  -32). 

To  obtain  Kelvin  (K)  readings,  use 

formula : 

K  =  (5/9)  (F  -32)  +  273.15. 
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SYMBOLS  AND  DEFINITIONS 


ai 
ar 
A,  B 
d 

dS 

8 

Hb 

Hi 

Ho 

Hr 

Hs 

Ht 

Hd 

Kr 

Kt 

k 

L 

Ho 

m 

n 

R 

Re 

T 

TP 

W 

a, a'  ,£ 

Y 

AJ. 

Prms 

e 


incident  wave  amplitude  at  a  spectral  line 
reflected  wave  amplitude  at  a  spectral  line 

real  and  imaginary  spectral  coefficients  from  an  FFT  analysis 
representative  armor  diameter  =  (W/y)1^3 
water  depth  at  the  toe  of  the  structure 
acceleration  due  to  gravity 

a  representative  breaking  wave  height  at  the  toe  of  the  structure 

incident  wave  height  (use  Hs  for  irregular  waves) 

deepwater  wave  height 

reflected  wave  height 

significant  wave  height 

transmitted  wave  height 

wave  dissipation  coefficient 

wave  reflection  coefficient 

wave  transmission  coefficient 

wave  number  =  2tt/L 

wavelength  at  the  toe  of  the  structure 

deepwater  wavelength  from  linear  theory  =  gT‘>/(27r) 

offshore  slope  seaward  of  the  structure 

number  of  layers  of  armor 

wave  runup 

Reynolds  number 

wave  period  (use  period  of  peak  energy  density  for  irregular  waves) 

period  of  peak  energy  density 

weight  of  armor  material 

empirical  wave  reflection  parameters 

specific  weight  of  armor  unit  material 

wave  gage  spacing 

average  root-mean-square  surface  water  level 
angle  of  the  seaward  structure  face 
kinematic  viscosity  of  water 
surf  similarity  parameter  =  tan  6// Hi/Lo 
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ESTIMATION  OF  WAVE  REFLECTION  AND  ENERGY  DISSIPATION 
COEFFICIENTS  FOR  BEACHES,  REVETMENTS,  AND  BREAKWATERS 


by 

William  IL  Seelig  and  John  P.  Ahrens 

I .  INTRODUCTION 

When  a  wave  encounters  a  coastal  structure  or  beach,  a  part  of  the  wave 
energy  is  dissipated.  The  remaining  energy  is  reflected  seaward  except  in  the 
case  of  a  permeable  or  overtopped  structure  (Fig.  1),  which  allows  transmission 
of  a  part  of  the  energy  to  the  leeward  side.  Wave  reflection  may  have  undesir¬ 
able  effects  because  the  reflected  waves  are  superimposed  on  the  incident  waves 
to  increase  the  magnitude  of  water  particle  velocities  and  water  level  fluctu¬ 
ations  seaward  of  the  structure.  These  enhanced  motions  may  be  a  hazard  to 
navigation  or  may  undesirably  alter  sediment  transport  patterns.  This  report 
presents  methods  for  estimating  wave  reflection  coefficients  for  beaches, 
revetments,  and  breakwaters  of  waves  approaching  the  structure  at  a  normal  angle 
of  incidence  (wave  crests  are  parallel  to  the  structure  axis). 

II.  LITERATURE  REVIEW 

Previous  investigators  have  experimentally  and  analytically  studied  wave 
energy  dissipation  and  reflection  characteristics  for  a  variety  of  structures. 
Various  prediction  techniques  have  been  proposed  to  estimate  reflection  coef¬ 
ficients  for  specific  types  of  energy  dissipation.  Miche  (1951)  proposed  a 
wave  reflection  coefficient  prediction  technique  that  is  often  quoted  in  lit¬ 
erature  (e.g.,  Sec.  2.54  in  L'.S.  Army,  Corps  of  Engineers,  Coastal  Engineering 
Research  Center,  1977).  He  assumed  that  there  is  some  critical  deepwater  wave 
steepness  below  which  the  reflection  coefficient  is  a  constant.  For  conditions 
where  wave  steepness  is  greater  than  the  critical  value,  the  reflection  coef¬ 
ficient  is  proportional  to  the  ratio  of  the  wave  steepness  to  the  critical 
value  of  wave  steepness.  Predictions  using  Miche’ s  approach  give  the  right 
order  of  magnitude  estimate  of  the  reflection  coefficient,  but  as  Ursell,  Dean, 
and  Yu  (1960)  illustrated,  predictions  may  be  conservative  by  a  factor  of  2. 

Moraes  (1970)  has  performed  some  of  the  most  extensive  laboratory  tests  to 
date  on  monochromatic  wave  reflection  from  a  variety  of  smooth  and  rough  slopes. 


- - 
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Madsen  and  White  (1976)  made  a  number  of  additional  carefully  controlled 
reflection  measurements  for  smooth  and  rough  steep-sloped  structures  under 
nonbreaking  wave  action.  Based  on  these  data,  they  developed  an  analytical- 
empirical  model  for  predicting  reflection  coefficients  for  rough  slopes  with 
nonbreaking  waves. 

Battjes  (1974)  used  Moraes'  data  to  develop  an  equation  for  predicting 
reflection  coefficients  for  smooth  slopes  where  the  slope  induces  wave  breaking 
This  technique  is  conservative  for  nonbreaking  (surging)  waves.  Ahrens  (1980) 
has  made  a  number  of  irregular  wave  reflection  coefficient  measurements  for 
overtopped  and  nonovertopped  plane  smooth  slopes. 

A  number  of  wave  reflection  measurements  for  laboratory  breakwaters  have 
been  made.  Seelig  (1980)  investigated  rubble-mound  and  caisson  breakwaters 
using  monochromatic  and  irregular  waves.  Brunn,  Gunbak,  and  Kjelstrup  (1979) 
measured  reflection  coefficients  for  rubble-mound  breakwaters  and  proposed  an 
empirical  prediction  technique.  Additional  breakwater  reflection  data  are 
available  in  Debok  and  Sollitt  (1978)  and  Sollitt  and  Cross  (1976).  Madsen 
and  White  (1976)  give  a  procedure  for  predicting  reflection  from  rubble-mound 
breakwaters  for  nonbreaking  waves. 

Chesnutt  and  Galvin  (1974)  and  Chesnutt  (1978)  have  made  some  of  the  most 
detailed  measurements  available  of  wave  reflection  from  laboratory  sand  beaches 
Little  prototype  data  are  available;  however,  Munk,  et  al .  (1963)  and  Suhayda 
(1974)  reported  reflection  measurements  for  beaches  exposed  to  extremely  low 
steepness  swell  waves. 


III.  EXPERIMENTAL  TECHNIQUES 


The  primary  emphasis  of  this  report  is  on  the  reanalysis  of  existing  data 
from  a  number  of  published  sources.  However,  some  additional  laboratory  data 
were  taken  to  supplement  the  sources;  these  data  are  reported  in  Appendix  A. 


Goda  and  Suzuki' s  (1976)  method  was  used  to  determine  wave  reflection  coef¬ 
ficients.  This  method  was  selected  because  with  the  test  setup  used  it  gave 
consistent  results  which  are  as  reliable  as  obtainable  with  other  currently 
used  procedures.  Experience  with  this  technique  suggests  that  the  error  is  on 
the  order  of  5  percent.  A  typical  wave  gage  setup  is  illustrated  in  Figure  2, 
and  a  detailed  discussion  of  the  analysis  method  given  in  Appendix  B.  The  test 
procedure  uses  three  gages,  located  a  minimum  of  6  meters  seaward  of  a  test 
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structure,  to  collect  simultaneous  wave  records  (incident  and  reflected  waves 
superimposed),  each  containing  4,096  data  points  at  a  sampling  interval  of  one- 
sixteenth  of  a  second.  A  fast  Fourier  transform  (FFT)  analysis  is  made  of  each 
record,  and  each  gage  pair  gives  an  estimate  of  the  reflection  coefficient  sub¬ 
ject  to  the  criteria  discussed  in  Appendix  B.  The  mean  of  the  three  estimates 
is  taken  as  representative  at  each  spectral  line,  and  an  energy-weighted  aver¬ 
age  is  taken  to  characterize  reflection  for  the  entire  spectrum  of  irregular 
waves.  The  significant  incident  wave  height,  Hs,  for  irregular  waves  (Goda 
and  Suzuki,  1976)  is  defined  as 

4  Urms 

Hs  ~  k  +  k£  ^ 

where  nr[ns  is  the  average  root-mean-square  (rms)  water  surface  displacement 
of  the  wave  records  at  the  three  gages,  and  Kr  the  reflection  coefficient. 

Data  collection  in  this  study  emphasized  obtaining  additional  data  on  wave 
reflection  on  smooth  slopes  and  examining  the  influence  of  one  or  more  layers 
of  armor  on  reducing  the  reflection  coefficient.  Monochromatic  and  irregular 
waves  were  tested. 

For  monochromatic  wave  conditions  (sinusoidal  wave  generator  blade  motion), 
the  wave  reflection  measurement  technique  was  slightly  modified.  The  wave¬ 
form  for  monochromatic  waves  is  described  by  a  Fourier  series  with  the  entire 
waveform  moving  at  the  speed  of  the  primary  wave  (Dr.  R.  Dean,  University  of 
Delaware,  personal  communication,  1980).  This  allows  the  wave  energy  appearing 
in  harmonics  of  the  primary  wave  to  be  considered  in  determining  the  reflec¬ 
tion  coefficient  (App .  B) . 

IV.  FACTORS  INFLUENCING  WAVE  REFLECTION 


The  conversion  of  wave  energy  concept  is  useful  for  defining  the  interre¬ 
lation  between  the  wave  reflection,  dissipation,  and  transmission  coefficients. 
Assuming  that  the  water  depth  remains  constant  seaward  and  leeward  of  the  struc¬ 
ture  the  partition  of  wave  energy  is  given  by 


1 


K£  + 


(2) 


where  Kr  is  the  reflection  coefficient,  the  ratio  of  wave  energy  lost 

through  dissipation  to  the  total  incident  wave  energy,  and  Kt  a  transmission 
coefficient  including  transmission  through  a  permeable  structure  and  trans¬ 
mission  by  overtopping  for  a  low-crested  structure.  In  an  idealized  monochro¬ 
matic  wave  situation  where  there  are  no  transfers  of  wave  energy  to  other  wave 
frequencies. 


and 


Ik 

Hi 


(3) 


Kt 


(4) 


where  H|,  Hr ,  and  Ht  are  the  incident,  reflected,  and  transmitted  wave  heights, 
respectively  (see  Fig.  1). 


II 


Rearranging  equation  (2)  gives 


Kr  =  V1  "  (Kd  +  Kf)  (5) 

which  clearly  shows  that  any  process  that  increases  the  sum  will  cause 

the  reflection  coefficient  to  decrease.  Figure  3  illustrates  equation  (5)  and 
the  nonlinear  relation  of  the  variables.  Note  that  for  a  given  value  of  the 
transmission  coefficient  the  reflection  coefficient  may  be  very  sensitive  to 
the  amount  of  energy  dissipation.  In  addition,  with  no  transmission  large 
values  of  energy  dissipation  will  allow  the  reflection  coefficient  to  be  rela¬ 
tively  large.  For  example,  with  90-percent  energy  dissipation  and  no  trans¬ 
mission,  the  reflection  coefficient  is  0.31  (see  Fig.  3). 


Pcf  Wove  Energy  Dissipated 


Figure  3.  Relation  between  wave  reflection,  transmission, 
and  dissipation  coefficients. 


V.  TYPES  OF  STRUCTURES  AND  RANGE  OF  CONDITIONS  TESTED 


Table  1  summarizes  the  sources  of  wave  reflection  coefficients  for  struc¬ 
tures  and  beaches  and  the  range  of  conditions  tested.  Three  types  of  structure 
are  considered:  smooth,  impermeable  slopes  with  no  overtopping;  revetments 
armored  with  one  or  more  layers  of  riprap  with  no  overtopping;  and  rubble-mound 
breakwaters  armored  with  stone  or  dolos. 


The  water  depth  at  the  toe  of  the  structure,  ds,  is  taken  as  a  character¬ 
istic  water  depth,  g  is  the  acceleration  due  to  gravity,  and  a  representative 
armor  unit  diameter,  d,  is  determined  from 


d 


(6) 


where  W  is  the  armor  weight,  and  y  the  specific  weight  of  the  armor  mate¬ 
rial.  A  measure  of  the  wave  breaker  height  that  could  occur  at  the  toe  of  the 
structure,  Hb,  is  given  by  Goda  (1975)  as 

Hb  =  0.17  Lpjl.O  -  exp  £-4 .712  ~  (l.O  +  15  mK  3  3  3jJ|  (7) 

where  L0  is  the  deepwater  wavelength  given  by  linear  wave  theory,  and  m  the 
tangent  of  the  slope  of  the  seabed  seaward  of  the  structure. 

Other  variables  summarized  in  Table  1  include  dimensionless  ratios  using 
Hj,  the  incident  wave  height  (significant  height  for  irregular  waves)  at  the 
toe  of  the  structure;  T,  the  wave  period  (period  of  peak  energy  density  for 
irregular  waves);  and  L,  the  wavelength  at  the  toe  of  the  structure. 


Only  those  tests  with  fully  turbulent  hydraulic  conditions  are  considered 
in  order  to  minimize  the  influence  of  viscous  effects  (Jonsson,  1966).  The 
Reynolds  number,  Re,  proposed  by  Madsen  and  White  (1976), 


_  R2  2n 
e  T  u  tan& 


(8) 


where  R  is  the  wave  runup  and  u  the  kinematic  viscosity  of  water  (about 
0.009  square  centimeter  per  second  at  20°  Celsius),  is  used  to  establish  which 
tests  are  fully  turbulent.  For  smooth  slopes  only  those  tests  with  Re  >  3  x  104 
are  analyzed;  for  rough  slopes  only  tests  with  Re  >  104  are  considered  (Jonsson, 
1966:  Madsen  and  White,  1976). 


VI.  TECHNIQUES  FOR  PREDICTING  REFLECTION  AND  ENERGY 
DISSIPATION  COEFFICIENTS 

Section  IV  showed  the  strong  dependence  of  the  magnitude  of  the  reflection 
coefficient  on  the  amount  of  wave  energy  dissipated  (also  on  the  amount  of  wave 
energy  transmitted  in  the  case  of  a  permeable  or  overtopped  structure) .  In  this 
section,  factors  that  influence  the  reflection  coefficient  are  systematically 
investigated,  and  empirical  prediction  formulas  are  developed.  Types  of  wave 
energy  dissipation  considered  include  losses  in  energy  due  to  structure-induced 
wave  breaking  and  wave  modification,  breaking  at  the  toe  of  a  structure  or 
in  the  surf  zone  seaward  of  the  structure,  structure  surface  roughness,  and 
internal  flow  in  permeable  sections  of  a  structure. 


1 


o  -o 

.=  c 

4-»  CQ  /->  sO 

Of  >£>  r- 

6  H  n 

UO>« 


<5  *3  -X  *H 

w  O  3  £ 

3  .c  n  3: 

U  W  3 

W  41  W  -C 

03  E  C 

(J  T3  B 

HI  C 

w  •*  a  a  c 

C  *3  O  C; 

Qj  O  «-4  <G  30 

■H  X  41  T)  ~ 

O  *■*  >  O  ts 

y  c  u  a: 

^  e  4i 

01  O)  T>  o  O 

o  a  4i 

O  O  -H  -3  -3 

•-1  u-  O  O 
C  <U  -rH  £  -C 

o  >  -a 

'ta  3  o  a*  n> 

*-*  O  E  E  E 


sine  blade  motion; 
irregular  waves. 


1 .  Modification  of  the  Wave  by  the  S t ructure  (Smooth  Slopes). 


For  a  structure  with  a  toe  water  depth-to-wave  height  ratio  greater  than 
five  and  wave  steepness  much  less  than  one-seventh,  the  interaction  of  the  wave 
and  structure  will  have  dominant  control  on  the  magnitude  of  the  reflection 
coefficient.  Miche  (1951)  proposed  that  the  reflection  coefficient  for  this 
situation  is  proportional  to  the  ratio  of  a  critical  wave  steepness  to  the  inci¬ 
dent  wave  steepness.  The  critical  steepness  is 


(£  In  (?) 


(9) 


where  Ho  is  the  deepwater  wave  height,  and  0  the  angle  the  structure  slope 
makes  with  the  horizontal,  in  radians.  Miche's  equation  gives  conservative 
results.  For  example,  it  overpredicts  monochromatic  wave  reflection  from  a  1 
on  15  slope  by  a  factor  of  2  (Ursell,  Dean,  and  Yu,  1960). 

Battjes  (1974)  recommends  the  equation, 


tan0 

;  ! 

»i 

\ 

Lo 

which  can  be  written  as 
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r 


0.1  tan' 0 


(10) 


(ID 


Battjes  (1974)  is  assuming  an  equation  similar  to  the  formula  proposed  by 
Miche  (1951)  where  the  critical  steepness  is 


(  =  0.1  tan20  (12) 

\  Lo/cr it 

This  criterion  gives  lower  and  more  realistic  values  of  the  reflection  coeffi¬ 
cient  than  Miche  (1951)  and  is  especially  useful  for  f,  <  2.3  where  breaking  is 
induced  by  the  structure  (for  plunging  breakers).  Figure  4  shows  the  compari¬ 
son  between  the  equations  of  Battjes  (1974)  and  Miche  (1951). 

The  following  revised  equation. 


Kr  =  tanh  (0.1  f/),  (13) 

is  recommended  to  give  a  conservative  prediction  of  reflection  coefficients. 

At  small  values  of  the  surf  similarity  parameter  (t,  ■  2.3), 

0.1  C?  ~  tanh  (0.1  f,2)  (14) 

and  equation  (13)  gives  the  same  results  as  equation  (10).  At  Jarger  values  of 
the  surf  similarity  parameter,  f, ,  equation  (13)  asymptotically  approaches  1.0 
and  gives  an  upper  bound  closer  to  the  data  than  equation  (10)  (see  Fig.  4). 


comparison  of  wave  reflection  coefficients  for  a  1  on  2.5  slope 
d  various  equations  to  predict  reflection  coefficients. 


An  improved  equation  for  predicting  reflection  coefficients  with  less  error 
in  the  estimates  is 


Kr  = 


£2  +  B 


1  + 


(15) 


where  a  and  B  are  empirical  coefficients  determined  from  the  laboratory 
data  (e.g..  Fig.  4).  The  value  of  B  increases  as  the  slope  becomes  flatter 
and  is  larger  for  irregular  waves  than  for  monochromatic  waves  (Fig.  5).  For 
slopes  with  cot0  <  6,  the  suggested  prediction  coefficients  are  a  =  1.0  and 


B  =  5.5  with  the  equation, 

Kr  = 

a  £2  | 

C2  +  B 

°r  1 

whichever 
is  smaller 

(16) 

K-j.  =  n  tanh  (0.1  £2) 
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Figure  5.  3  as  a  function  of  structure  slope. 

2 .  Breaking  at  the  Toe  or  Seaward  of  the  Structure. 

If  the  water  depth  at  the  toe  of  the  structure  is  less  than  five  times  the 
incident  wave  height  or  if  the  wave  steepness  is  large,  significant  additional 
wave  energy  loss  may  result  from  wave  steepness/water  depth-limited  breaking. 
The  dimensionless  ratio  describing  this  type  loss  is  the  ratio  of  the  incident 
wave  height  to  the  maximum  possible  breaker  height,  (H^/Ht,) ,  where  Hj,  is 
given  by  equation  (7).  This  ratio  includes  the  influence  of  offshore  slope, 
water  depth  at  the  toe  of  the  structure,  and  wave  steepness,  and  gives  a  meas¬ 
ure  of  breaking  at  the  toe.  The  suggested  empirical  coefficient  to  account 
for  this  type  energy  loss  in  predicting  reflection  coefficients  is 

*  ■  “»[-  °-5(|)13] 

for  use  with  equation  (16),  where  a  is  a  reflection  coefficient  reduction 
factor . 
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•  Irregulor  Waves  (Ahrens,  1980) 

O  Monochromatic  Woves  (  Ursell ,  Deon ,  and  Yu, 

I960;  Moraes,  1970;  This  Study  ) 
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(17) 


3.  Influence  of  Surface  Roughness. 

Armor  units  placed  on  the  surface  of  a  smooth  structure  will  increase  the 
amount  of  energy  loss  in  a  wave  encountering  the  structure,  thereby  reducing 
the  amount  of  wave  reflection.  The  suggested  prediction  equation  for  a  revet¬ 
ment  with  one  layer  or  armor  rock  with  representative  diameter,  d,  is 


for  use  with  equation  (16),  where  L  is  the  wavelength  at  the  toe  of  the  struc¬ 
ture.  This  equation  was  developed  from  the  data  in  Table  1. 

Figure  6  illustrates  the  joint  influence  of  a  relative  armor  roughness 
parameter,  /d/L  cotG,  and  a  relative  breaking  height  parameter,  on  the 

reflection  coefficient  reduction  factor,  a.  An  examination  of  equation  (18) 
and  Figure  6  indicates  that  if  all  other  factors  remain  fixed,  the  reflection 
coefficient  will  decrease  as  the  ratio  of  the  stone  size  to  wavelength,  d/L, 
increases,  as  the  cot0  increases  (the  slope  becomes  flatter),  or  as  the  ratio 
of  the  incident  wave  height  to  the  breaking  wave  height, (H^/H^) ,  increases. 
Figure  7  shows  a  comparison  between  predicted  reflection  coefficients  using 
equations  (18)  and  (16)  versus  observed  reflection  coefficients  for  monochro¬ 
matic  and  irregular  waves  on  a  1  on  2.5  slope  armored  with  one  layer  of  stone 
with  d/ds  =  0.15.  The  correlation  coefficient  is  0.98  for  monochromatic  waves 
and  0.94  for  irregular  waves. 


The  ratio  of  armor  stone  diameter  to  incident  wave  height,  d/H^ ,  on  the 
average  has  little  influence  on  the  reflection  coefficient  for  one  layer  of 
armor,  so  this  parameter  is  not  included  in  equation  (18).  Some  deviation 
from  equation  (18)  occurs  where  stone  size  is  much  larger  than  wave  height  and 
resulting  predictions  are  conservative.  For  example,  where  the  stone  size-to- 
wave  height  ratio  is  greater  than  2.0,  equations  (16)  and  (18)  overpredict 
reflection  coefficients  by  an  average  of  6  percent. 


./a/L  cot 9 

Figure  6.  Joint  effect  on  one  layer  of  armor  and  H^/H^  on  the 
reflection  coefficient  reduction  factor,  a. 
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4.  Influence  of  Multiple  Layers  of  Armor. 


As  the  number  of  layers,  n,  of  armor  on  a  revetment  increases,  the  amount 
of  wave  energy  dissipated  increases  and  the  reflection  coefficient  decreases. 

In  addition,  as  the  size  of  the  stone  increases  relative  to  the  wave  height, 
the  roughness  becomes  more  effective  and  the  reflection  coefficient  decreases. 
Table  2  gives  selected  Values  of  a  correction  factor,  a',  where 


a 


exp 


r  cote 


(19) 


Table  2.  Correction  factor  due  to  multiple 
_ layers  of  armor.1 _ 

a '  _ 


L  n 

d/Ht 

Two 

Three 

Four 

<0.75 

0.93 

0 .88 

0.78 

0.75  to  2.0 

0.71 

0.70 

0.69 

>2.0 _ _ 

0.58 

...  °.-52.J 

0.49 

'cote  =  2.5,  d/ds  =  0.15,  0.004  <  ds/gT2 
<  0.03. 


for  multiple  layers  of  armor.  These  coefficients  were  obtained  by  taking  the 
average  of  the  ratios  of  the  measured  reflection  coefficients  for  two,  three, 
and  four  layers  of  armor  to  predicted  coefficients  for  a  slope  with  one  layer 
of  armor.  Only  one  slope,  cot0  =  2.5,  and  stone  size-to-water  depth  ratio, 
d/ds  =  0.15,  was  tested. 

5 .  Wave  Reflection  from  Sand  Beaches. 

Chesnutt  (1978)  has  the  most  extensive  data  set  of  wave  reflection  coeffi¬ 
cients  from  laboratory  sand  beaches.  Unfortunately,  there  are  little  prototype 
data  available.  Chesnutt  and  Galvin  (1974)  and  Chesnutt  (1978)  found  that  many 
factors  influence  the  magnitude  of  the  reflection  coefficient.  Their  data 
suggest  that 


Kr  =  ~~ —  ;  B  =  5.5  (20) 

S2  +  B 

can  be  used  to  estimate  reflection  coefficients  with  the  beach  slope  at  the 
Stillwater  level  intercept  used  to  determine  £.  Use  a  =  1.0  for  conservative 
estimates  of  Kr  and  a  =  0.5  to  give  predictions  of  the  average  reflection 
coefficient  measured  throughout  a  test  (Fig.  8). 


Kr  0.3 


Mean  Kr 


Range  of  Kr 


Figure  8.  Wave  reflection  coefficients  from  laboratory 
beaches  (from  Chesnutt,  1978). 

6.  Rubble-Mound  Breakwaters. 

An  upper  limit  or  conservative  estimate  of  Kr  for  breakwaters  armored 
with  rock  or  dolos  may  be  obtained  using 

Kr  =  vHLa  ’>  a  =  °’6«  3  "  6-6  (21) 

K  +  p 

Ninety-five  percent  of  all  observed  laboratory  breakwater  wave  reflection 
coefficients  fall  below  this  prediction  equation  for  data  sets  c,  d,  g,  and  h 
outlined  in  Table  1. 

More  reliable  predictions  of  wave  reflection  coefficients  for  rubble-mound 
breakwaters  may  be  made  using  the  method  of  Madsen  and  White  (1976)  (also  see 
Seelig,  1979).  Equations  (16)  and  (18)  should  be  used  with  the  Madsen  and 
White  (1976)  method  to  estimate  energy  dissipation  on  the  seaward  face  of  the 
breakwater  caused  by  the  outer  layer  of  armor  units.  Figure  9  shows  sample 
laboratory  measurements  (Sollitt  and  Cross,  1976)  and  predicted  reflection  and 
transmission  coefficients  for  a  rubble-mound  breakwater.  Observed  and  predicted 
reflection  coefficients  have  the  best  agreement  for  wave  conditions  in  the  tur¬ 
bulent  zone,  but  deviate  where  the  Reynolds  number  becomes  less  than  10^  due 
to  laboratory  scale  effects. 


*’>*>»■  '  • 


Figure  9.  Predicted  rubble-mound  breakwater  wave  reflection 
and  transmission  coefficients  (laboratory  data 
fror  Sollitt  and  Cross,  1976), 
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7 .  Spectral  Resolution  of  Wave  Reflection. 


The  significant  wave  height  and  period  of  peak  energy  density  are  used  to 
characterize  irregular  wave  conditions  in  this  report.  However,  a  more  detailed 
analysis  shows  that  the  reflection  coefficient  varies  as  a  function  of  wave 
frequency  for  irregular  waves.  Figure  10  illustrates  the  decrease  in  reflection 
coefficient  as  a  function  of  wave  frequency  that  is  typical  of  waves  breaking 
on  a  smooth  impermeable  1/2  slope  (£  <  2.3).  Nonbreaking  waves  have  a  different 
characteristic  shape  of  the  reflection  coefficient  as  a  function  of  wave  fre¬ 


quency.  Kr  increases  as  a  function  of  f  for  frequencies  higher  than  the 


frequency  of  peak  energy  density  (Fig.  11).  The  shift  to  high  frequencies  seems 
to  occur  because  wave  energy  is  transferred  from  low  to  higher  frequencies  due 
to  nonlinear  effects  when  the  waves  interact  with  the  structure.  Note  that  this 
energy  shift  may  produce  a  range  of  wave  frequencies  in  which  more  wave  energy 
is  moving  away  from  the  structure  than  is  incident  to  the  structure,  and  the 
local  reflection  coefficient  may  be  larger  than  1.0  over  this  range  of  fre¬ 
quencies.  Caution  should  be  used  when  trying  to  obtain  information  from  the 
highest  frequency  part  of  the  spectrum  above  approximately  the  95-percent  cumu¬ 
lative  energy  density  level  because  the  signal-to-noise  ratio  is  low  and  the 
wave  speed  is  poorly  known  (Mansard  and  Funke,  1979). 


8 .  Reflection  Coefficient  Prediction  Equations. 


Table  3  summarizes  the  equations  recommended  for  estimating  reflection 
coefficients  for  slopes,  revetments,  rubble-mound  breakwaters,  and  beaches. 
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Figure  10.  Wave  reflection  coefficient  as  a  function  of 

wave  frequency  for  an  irregular  wave  condition 
with  breaking  waves. 
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Figure  11.  Wave  reflection  coefficient  as  a  function  of  wave 

frequency  for  an  irregular  nonbreaking  wave  condition 


VII.  EXAMPLE  PROBLEMS 


The  following  example  problems  illustrate  the  methods  of  predicting  reflec¬ 
tion  coefficients  presented  in  this  report. 


**************  * EXAMPLE  PROBLEM  1*  *************** 

GIVEN:  A  smooth  impermeable  revetment  (nonover topped)  has  a  toe  water  depth, 
ds  =  7.62  meters,  a  slope  cot0  =  2.0,  and  the  offshore  slope  is  m  =  0.02. 

FIND:  The  wave  reflection  coefficient  and  fraction  of  wave  energy  dissipated 
for  a  wave  with  =  3.05  meters  and  T  =  10  seconds.  Illustrate  the  influ¬ 
ence  of  wave  height  and  period  on  1C,,  and  show  the  effect  of  reducing  the 
slope  to  cot6  =  5. 

SOLUTION:  From  equation  (7), 


Hb  =  0.17  L0  j  1.0  -  exp  [-  4.712  ^  ^1  +  15  m  1-33)]) 

Hb  =  0.17  (  1.56  x  102)jl  -  exp  [-  4.712  ^  (l  +  15  (0 . 02) 1 • 3 3)]} 
From  equation  (17) 

[-°'5(sf)  J’  ]-  0.807 


5.85  m 


a  =  exp 


From  equation  (10) 


tanO 


0.5 

3.05 

156 


3.58 


and  from  equation  (15) 


a  t,2  =  0.807  (3. 58)2 

f.2  +  3  (3 . 58)  2  +  5.5 


0.56 


The  energy  dissipation  coefficient  for  this  example  is  =  0.69,  or  69  per¬ 
cent  of  the  incident  wave  energy  is  dissipated  (from  Fig.  3).  Other  reflection 
coefficient  calculations  for  5-,  10-,  and  20-second  periods  for  wave  heights 
between  0.3  and  4.4  meters  are  summarized  in  Figure  12.  Predictions  are  also 
shown  for  a  structure  with  cote  =  5.  Figure  12  illustrates  the  influence  of 
wave  height,  period,  and  structure  slope  on  Kr. 
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Figure  12.  Predicted  wave  reflection  coefficients  for  smooth 
impermeable  slopes  with  no  overtopping. 

**************  * EXAMPLE  PROBLEM  2**************** 

GIVEN :  The  wave  conditions  in  example  problem  1. 

FIND:  The  wave  reflection  coefficients  if  one  layer  (n  =  1)  or  two  layers 

(n  =  2)  of  4,500-kilogram  (5  tons)  rock  at  2,700  kilograms  per  cubic  meter 
(169  pounds  per  cubic  foot)  were  added  as  armor  to  the  revetment  with  cotB  = 
2.0. 

SOLUTION :  The  armor  material  in  this  example  has 


i  /  i  l  /  i 

'(?)  -M 


using  equation  (6).  For  the  case  of  T  =  10  seconds  and  H  =  3.05  meters 
equation  (18)  gives 


a  =  exp 


[-1. 7$  cote  -  0.5  (|l*j  ] 


=  0.536 


and  from  equation  (16) 


=  a  ^  =  0. 536(3. 58); 

f  +  f.2  ~  5.5  +  (3 . 58)  ■ 


=  0.37 


The  energy  dissipation  coefficient  from  Figure  2  is  Kj  =  0.86,  86-percent  dissi¬ 
pation  or  17  percent  more  dissipation  than  occurred  for  the  smooth  slope  (see 
example  problem  1).  Sample  predicted  reflection  coefficients  are  given  in 
Figure  13.  The  preliminary  information  in  Table  2  suggests  that  further  re¬ 
duction  in  the  reflection  coefficients  could  be  achieved  by  adding  a  second 
layer  of  armor  (n  =  2)  for  wave  heights  less  than  3  meters  (Fig.  13). 


Hj  (m) 

Figure  13.  Wave  reflection  coefficients  for  a  smooth  revetment 

and  revetments  with  one  and  two  layers  of  armor  stone. 


VIII.  SUMMARY 

Methods  for  predicting  wave  reflection  and  dissipation  coefficients  for 
beaches,  nonover topped  revetments,  and  breakwaters  are  presented.  Types  of 
wave  energy  dissipation  considered  are  wave  breaking  induced  by  the  structure, 
wave  breaking  at  the  toe  of  the  structure,  turbulence  produced  by  wave  inter¬ 
action  with  the  outer  layer  of  armor,  and  flow  through  additional  layers  of 
armor.  These  techniques  are  combined  with  the  method  of  Madsen  and  White 
(1976)  to  estimate  reflection  and  transmission  coefficients  for  permeable 
rubble-mound  breakwaters.  Factors  considered  when  making  reflection  coeffi¬ 
cient  estimates  include  structure  slope,  water  depth  at  the  toe  of  the  struc¬ 
ture,  offshore  slope,  incident  wave  height  and  period,  the  size  and  number  of 
layers  of  armor  units,  and  the  type  of  structure.  Techniques  presented  apply 
to  breaking  and  nonbreaking  (surging)  waves  and  can  be  used  for  monochromatic 
and  irregular  wave  conditions. 
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APPENDIX  A 


LABORATORY  WAVE  REFLECTION  DATA 

This  appendix  includes  tables  of  wave  reflection  data  (Tables  A-l  to  A-7) 
obtained  as  a  part  of  this  study.  The  following  variables  are  used: 

ID  -  an  identification  code  assigned  to  each  data  run 

H  -  the  incident  wave  height  (centimeter);  the  significant  wave  height 
for  irregular  waves 

T  -  the  wave  period  (second),  the  period  of  peak  energy  density  for 
irregular  waves 

SURF  -  the  surf  similarity  parameter  =  tan6//H^/gT2 

H/HB  -  the  incident  wave  height  divided  by  the  maximum  breaker  height 
expected  at  the  toe  of  the  structure 

D/H  -  water  depth  divided  by  incident  wave  height 

KR  -  reflection  coefficient 

QP  -  the  spectral  peakedness  parameter  for  irregular  wave  conditions 
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OTCBttlD  PAOK  BUNK-MOT  FILMED 


Tab 1 e  A-l .  Wave  reflection  from  a  1/15.0  smooth  slope  (monochromatic  waves)_. 
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Table  A-2 .  have  reflection  from  a  1/2.5  smooth  slope  (monochromatic  waves)^ 
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Table  A-3.  Wave  reflection  from  a  1/2.5  slope  with  one  layer  of  armor 
_ (monochromatic  waves). _ 
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Table  A-4 .  Wave  reflection  from  a  1/2.5  slope  with  two  layers  of  armor 
_ (monochromatic  waves). _ _ 
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Table  A-5.  Wave  reflection  from  a  1/2.5  slope  with  three  layers  of  armor 
(monochromatic  waves). _ _ _ _ _ — — 
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Table  A-6 .  Wave  reflection  from  a  1/2.5  slope  with  four  layers  of  armor 
(monochromatic  waves)  . _  _ _ — - 
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APPENDIX  B 


METHOD  OF  MEASURING  WAVE  REFLECTION  COEFFICIENTS 

The  method  of  Goda  and  Suzuki  (1976)  is  used  to  determine  laboratory  reflec¬ 
tion  coefficients  for  monochromatic  and  irregular  conditions.  Also  used  is  the 
energy  balance  approach  for  both  types  of  waves,  so  that  wave  energy  transfer 
between  frequencies  and  variable  amounts  of  reflection  over  a  range  of  fre¬ 
quencies  can  be  considered.  This  approach  gives  a  reflection  coefficient  that 
is  formally  defined  as  the  square  root  of  the  ratio  of  the  reflected  wave  energy 
to  incident  wave  energy.  For  an  idealized  case  where  no  energy  transfers  occur, 
the  reflection  coefficient  is  the  ratio  of  reflected  and  incident  wave  heights. 
Reflection  coefficients  are  determined  by  placing  two  or  more  gages  several 
wavelengths  seaward  of  the  structure.  Each  pair  of  gages  then  gives  an  esti¬ 
mate  of  reflection  coefficients. 

In  these  experiments  wave  records  were  sampled  simultaneously  at  three  wave 
gages  (Fig.  2)  at  a  rate  of  16  times  a  second  to  obtain  4,096  data  points  for 
each  run.  An  FFT  was  then  performed  on  each  wave  gage  record  to  determine  real 
and  imaginary  spectral  coefficients,  A  and  B,  at  each  spectral  line  j.  Let  the 
subscripts  j  and  2  indicate  the  landward  and  seaward  gages  in  a  pair.  The  re¬ 
flected  and  incident  wave  amplitudes  for  each  gage  pair  for  each  spectral  line 
are  then  given  by 

1  j A2  -  A[  cos  k.A£  -  Bj  sin  k.Af.)2  +  (B2  +  A)  sin  k A£  -  B]  cos  kAt)2  (B-l) 

2l  sin  kAtl 


1  A A  2  -  A;  cos  k&i  +  Bj  sin  tcAl)2  +  (Bj  -  Aj  sin  kfil  -  Bi  cos  k&i.)2  (B— 2) 

ar  '  YTiTTW 


A,B  =  spectral  coefficients 

k  =  wave  number  =  —•  (B-3) 

M  =  gage  spacing 

Only  gage  pairs  with 

0.05  £  ~  <  0.45  (B-4) 

are  used  in  the  analysis,  and  wavelength,  L,  is  determined  from  linear  theory 
for  irregular  waves, 

L  ■ ta"h  (¥)>  «-5> 

and  may  be  found  using  Dean's  (1974)  stream-function  theory  for  steep  monochro¬ 
matic  waves  (see  App.  C) . 

All  estimates  of  reflection  coefficients  found  using  the  above  procedure 
are  averaged  at  each  spectral  line  to  give  an  incident  wave  amplitude  and  re¬ 
flection  coefficient  for  line  j: 
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=  average  incident  wave  amplitude  at  line  j 


(ai)j 

(Kr)j 


average  reflection  coefficient  at  line  j 


The  reflection  coefficient  is  then  determined  by  taking 


(B-6) 


Irregular  wave  information  is  displayed  in  the  form  of  band  spectra,  using  11 
lines  per  band  and  using  a  variation  of  equation  (B-6)  to  determine  the  reflec¬ 
tion  coefficient  for  each  band. 

In  the  case  of  monochromatic  waves,  a  nonlinear  waveform  is  described  by 
a  Fourier  series  with  each  component  moving  at  the  speed  of  the  primary  wave, 
and  equation  (B-6)  is  used  to  determine  the  reflection  coefficient. 


APPENDIX  C 


NONLINEAR  WAVELENGHTS  AND  WAVE  SPEED 

In  the  real-time  analysis  of  wave  reflection  it  is  necessary  to  know  the 
wavelength  or  wave  speed.  Linear  theory  gives  excellent  predicitons  for  low 
steepness  waves,  but  tends  to  underestimate  both  length  and  speed  for  large 
waves . 


Dean  (1974)  gives  tabular  values  of  wave  speed  and  wavelength  for  finite 
height  waves  that  can  be  approximated  by  the  empirical  relation, 


where  L  and  C  are  wave  speed  and  wavelength,  LQ  and  C0 
speed  and  wavelength  determined  from  linear  theory  where 

L  -  gT2 
L°  "  2  IT 


are  deepwater  wave 


(C-2) 


L^  is  the  local  length  determined  from  linear  or  Airy  theory  and  a  and  b 
are  empirical  coefficients.  Airy  wave  theory  predictions  and  values  of  a  and 
b  are  plotted  as  a  function  of  ds/L0  in  Figure  C-l,  where  dg  is  the  water 
depth. 


Figure  C-l.  Coefficients  for  approximating  nonlinear  wave  speed 

and  wavelength  determined  from  stream-function  theory. 
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